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ABSTRACT 

The M31 X-ray source RX J0042. 3+4115 was originally identified as a black hole binary 
because it displayed characteristic low state variability at conspicuously high luminosities; un- 
fortunately, this variability was later found to be artificial. However, analysis of 84 Chandra 
ACIS observations, an HST ACS/WFC observation, and a 60 ks XMM- Newton observation has 
supplied new evidence that RXJ0042. 3-1-4115 is indeed a black hole binary. The brightest optical 
star within 3ct of the position of RXJ0042. 3+4115 had a F435W (-B) magnitude of 25.4+0.2; 
Mb > —0.4, hence we find a low mass donor likely. RX J0042. 3+4115 was persistently bright 
over ~12 years. Spectral fits revealed characteristic black hole binary states: a low/hard state at 
2.08+0.08 xlO^** erg s"\ and a steep power law state at 2.41+0.05x10^® erg s~^ (0.3-10 keV). 
The high luminosity low state suggests a ~20 Mq primary; this is high, but within the range of 
known stellar black hole masses. The inner disk temperature during the steep power law state 
is 2.24+0.15 keV, high but strikingly similar to that of GRS 1915+105, the only known Galactic 
black hole binary with a low mass donor to be persistently bright. Therefore RX J0042. 3+4115 
may be an analog for GRS 1915+105; however, other mechanisms may account for its behavior. 
We find compelling evidence for an extended corona during the steep power law state, because 
compact corona models where the seed photons for Comptonization are tied to the inner disc 
temperature are rejected. 

Subject headings: x-rays: general — x-rays: binaries — black hole physics 



Introduction 



The M31 X-ray sour c e RX J0042.3+4115, named 
following Supper et al. ( 1997[ ). was o riginally iden 



X-ray binaries, whether the accret or is a neu- 
tron star or black hole (see e.g. | van der Khi 



tified as a black hole X-ray binary bv lBarnard et al 
(|2003l ). after analyzing 4 XMM- Newton obser- 
vations from 2000-2002. It apparently exhib- 
ited power density spectra (PDS) that were 
well described by a broken power law, with 
spectral index a changing from ~0 to ^1 at 
some break frequency; such a PDS is consis- 
tent with the low/hard states observed in all 



1994 IWiinands fc van der Klii Il999r Neutron 
star X-ray binaries tend to exhibit such be- 
haviour at luminosities around lO'^^-lO'^'' erg 
s"\ yet RXJ0042. 3+4115 exhibited this vari- 
ability at 0.3-10 keV luminositi es of ^1-3x10^® 

(|2003l ) concluded that 



Barnard et al. 



erg s 

RX J0042.3+4115 contained a black hole. 

However, it was later discovered that these 
PDS and those reported by other groups were 
contaminated by artifacts caused by the XMM- 
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Newt on data reduction software ([Barnard et al 
20071 ). The problem arose because all XMM- 



newton liglitcurves start at the arrival time of 
the first photon by default; hence, source and 
background lightcurves, and lightcurves from the 
three EPIC detectors— MOSl, M0S2, and pn— 
are asynchronous by default. Combining these 
lightcurves (e.g. combining instruments, or back- 
ground subtraction) often resulted in PDS with 
artificial broken power law shapes. 

Therefore we make no use of the PDS and in- 
stead rely on our well established method of us- 
ing low state emission spectra (power law emission 
with photon ind ex 1.4-1.7 and little to no thermal 
emission, iMcClintock fc Rcmillard l2006l) at con- 
spicuously high luminosities to identify black hole 
candidates ('Barnard et al.l l2008t iBarnard fc Kolbl 
[2009: Barnard et al. .2011). We present the most 
detailed iustification of our selection criteria in 



Barnard et all ( 20111) 



In this paper we present our analysis of 84 
Chandra ACIS observations of RXJ0042. 3-1-4115 
over ~12 years, and our serendipitous HST obser- 
vation, along with our re-analysis of the 60 ks 2002 
XMM-Newton observation. We use the HST data 
to place RX J0042.3+4115 in M31, and argue for a 
low mass donor. We use long-term and short-term 
variability, and also emission spectra, to reinstate 
RXJ0042.3-H4115 as a black hole candidate. We 
discuss the observations and data analysis in the 
next section, followed by our results in Section 3, 
and a discussion in Section 4. 

2. Observations and analysis 

We analyzed 84 Chandra ACIS observations of 
the central region of M31, spaced over '^12 years, 
using CI AO version 4.3. For each observation 
we extracted 0.3-7.0 keV source and background 
spectra from circular regions with 10" radius; the 
background region was close to the source region, 
and source free. Corresponding response matrices 
and ancillary response files were also made. We 
obtained 0.3-10 keV luminosities from each obser- 
vation using XSPEC version 12.6.0. 

Observations with >200 net source counts were 
freely fitted with absorbed power law models; 
spectra were grouped to give at least 20 counts 
per bin. For observations with <200 net source 
counts we assumed an absorbed power law model 



with iVn =1.0x10^^ atom cm~^ and F = 1.5, and 
found the 0.3-10 keV luminosity equivalent to 1 
count s^^, then multiplied this conversion factor 
by the intensity; we chose this model because it 
approximates the best fit to our deepest Chandra 
observation of RXJ0042. 3+4115 in its low state. 
Luminosity uncertainties for freely fitted spectra 
are estimated by XSPEC by calculating a range 
of fluxes obtained by varying the emission param- 
eters; the uncertainties for the faint spectra are 
derived directly from intensity uncertainties. 

Additionally, we analyzed the 60 ks 2002 Jan- 
uary XMM-Newton observation of M31 (Rev 381) 
with SAS version 10.0.0. We extracted 0.3-10 
keV EPIC-pn lightcurves and spectra from circu- 
lar source and background regions with 15" radius. 
The background region was on the same CCD as 
the source region, and at a similar ofF-axis angle. 
The spectra were grouped to ensure a minimum of 
50 counts per bin. 

RXJ0042. 3-1-4115 was serendipitously observed 
in one of our HST ACS/WFC observations of M31 
transients. Observation j9udl7010 was made on 
2009 August 25, with the F435W fiher for 4360 
s. We registered a combined Chandra 0.3-7.0 
keV image and the HST image to the same B 
band i mage provided by th e Local Group Galaxy 
Survey ([Massev et al.[ [20061 ) using the IRAF task 
CCMAP; we used X-ray bright globular clusters to 
register the Chandra image, and bright, unsatu- 
rated stars to register the HST image. We deter- 
mined the best X-ray position of RX J0042. 3-1-4115 
with the IRAF task imcentroid. Wc used the 
IRAF package DAOPHOTII to perform photom- 
etry on the stars within 3a of the position of 
RXJ0042.3+4115. 

3. Results 

3.1. The search for an optical counterpart 

The centroid of the X-ray emission from 
RX J0042.3-h4115 was located at RA = 00:42:22.954, 
Dec = 41:15:35.23, with la uncertainties of 0.009" 
in RA and 0.007" in Dec. Combining this with 
the r.m.s. uncertainties in registration yields X- 
ray positional uncertainties of 0.09" in RA and 
0.19" in Dec. 

Figure [Tj shows a detail our HST image, su- 
perposed with an ellipse representing the 3a un- 
certainties in the position of RX J0042.3-f4115. 
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There are several stars within the ehipse, the 
brightest of which has a Vega B magnitude of 
25.4±0.2. We therefore constrain the B band mag- 
nitude to >24.8. We see no evidence for a back- 
ground galaxy; the region is relatively uncrowded, 
as it is from the M31 bulge, and the de- 
tection limit is B < 28. There is no evidence 
for a counterpart in th e 2MASS All Sky Catalog 
(|Skrutskie et all l2006[ ). hence RXJ0042. 3+4115 
cannot be associated with a late type star in our 
Galaxy. We conclude that it is located in M31. 

The distance modulus for M31 '^24.5, and we 
can estimate Ab by using the e mpirical relations 
A^iVn/l-SxlO^i atom cm -2 (jPredehl fc SchmittI 
19951 ). and E(B-V) - Av/3. The column density 



varied significantly between observations, hence 
the variable component was probably internal to 
the system. Since the donor is unlikely to suffer 
this extra absorption, we assume A^h — 1.0x10^^ 
atom cm~^ (see below), and = 0.7. Therefore 
Mb > -0.4. 

The known counterparts of high mass X-ray 
binaries (HMXBs) in the SMC have apparent V 
magnitudes in the range 13 < my < 18, and 
B -V in the range -0.32 < B-V < .06 (see e.g. 
Coe et al.|[200l lAntoniou et"aLll2009[) . For a dis- 
tance of '^60 kpc, this equates to —6 < Mb ^ — 1, 
all brighter than our threshold of Mb > —0.4. 

The three known BH HMXBs are Cygnus X-1, 
LMC X-1 and LMC X-3. Cygnus X-1 has a coun- 
terpart withMv^- 6. 5 (fWalborn 1973 ) and B-V 
= 0.8 ([HiltneilEgHel ). The counterpart to LMC 
X-1 has m agnitude V = 14.6 0ib0.02 and B-V = 

0.17±0.08 (jOrosz et al.ll2009h : hence Mb 4.2 

for a distance of 50 kpc. LMC X-3 has a B '-^ 17 
counterpart (see e.g. iBrocksopp et al.l 120011 ) . and 
Mb ~ —1.5; this is ~8cr brighter than the bright- 
est star within the ellipse. 

RXJ0042. 3-1-4115 is therefore most likely to be 
a low mass X-ray binary (LMXB). 

3.2. Time variability 

We present the ^12 year 0.3-10 keV luminosity 
lightcurve of RX J0042.3-b4115 created from the 
84 ACIS observations in Fig. [2 circles mark lu- 
minosities from freely-fitted bright spectra, while 
crosses represent faint spectra where we assume 
A^H = 1.0x10^1 atom cm s"^, and F = 1.5. Un- 
certainties are quoted at a la level. The lightcurve 



is extremely variable, with the luminosity varying 
over ~0. 5-3x10^^ erg s~^. We note that the spec- 
tral fits to all bright observations are consistent 
with a constant F, although Nh varied by a factor 
^^5; the faint observations may be up to 40% times 
brighter. 

RX J0042. 3+4115 appears to be persistently 
bright. By contrast, most Galactic black hole 
LMXBS are transient. One exception is GRS 
1915+105; a 7 year RXTE/ASM lightcurve of 
GRS 1915+1 05 showed it to be persistently bright 
(jMcClintock fc Remillard .2006). Another possi- 
ble exception is GRS 1758-258; it is thought to be 
a Galactic LMXB, but its true nature is not con- 
firm ed due to the high degree o f absorption (see 
e.g. iMuiioz-Arionilla et al.l boiol and references 
within) 

We also examined the short term variability of 
RXJ0042.3+4115 during the long XMM-Newton 
observation. We present the 0.3-10 keV EPIC-pn 
intensity lightcurve for RX J0042.3+4115 in Fig.H 
along with the background lightcurve in grey for 
comparison. The intensity varies by a factor 2 
(4(7 deviation) in <10 ks, hence the emission is 
dominated by a single source. The probability 
of RXJ0042. 3+4115 consisting of muhiple bright 
variable sources is very low, especially since it is 
not associated with any globular cluster. 

3.3. Spectral analysis 

3.3.1. Chandra observation OBSID1575 

The longest ACIS observation of RX J0042.3+4115 
was OBSID1575, an ACIS-S observation with a 
■~40 ks exposure time; the net source spectrum 
contained 7690 photons. An on-axis source with 
this intensity would be in danger of pile-up; how- 
ever, RX J0042.3+4115 was -4' off-axis, and the 
photons were spread over a large number of pixels 
(>100). Each incoming photon is assessed by its 
impact on a 3x3 array of ACIS pixels ; "good" 
photons are detected in only 2 of the 9 pixels, 
while cosmic rays etc. are detected in 3 or more 
(|Davisll200l[ ). We therefore estimated the proba- 
bility of pile up from the brightest pair of pixels; 
this pair accumulated 807 photons over ^^38 ks, or 
one photon every ^14 frames. Hence, we conclude 
that pile up is unlikely to have been significant. 

The 0.3-7.0 keV spectrum of RXJ0042. 3+4115 
during observation OBSID1575 is well described 



3 



by an absorbed power law, with line-of-sight ab- 
sorption A'^H = 1.0±0.2xl0^^ atom cm~^ and pho- 
ton index r 1.46±0.05; xVdof = 211/204. The 
0.3-10 keV luminosity was 2.08±0.08x 10^^ erg 
s~^. Uncertainties are quoted at a 90% confidence 
level. Figure |4] shows the unfolded 0.3-7.0 keV 
spectrum multiplied by the channel energy, assum- 
ing the best fit absorbed power law model. 

When a disk blackbody component was added 
to the power law emission, XSPEC set the inner 
disk temperature to 8.2xl0~^ keV, with A'h ~ 
1.0±0.2xl02i atom cm-\ and T = 1.46±0.05; 
X^/d.o.f = 211/202. Hence there is no trace of 
a disk component in the 0.3-7.0 keV spectrum. 
We therefore conclude that RX J0042. 3-1-4115 was 
in its low state during this observation; since the 
theoretical upper luminosity limit for low states 
in neutron star X-ray binaries is ~3xl0^^ erg s^^ 
RXJ0042.3 +4115 is a likely bl ack hole candi- 

and references 



folded spectrum multiplied by channel energy, 
assuming the best fit model. The 0.3-10 keV 
luminosity was 2.41±0.05xl0^* erg s~^, with 
the power law component contributing ~45%. 
Such a spectrum is consistent with the steep 
power law black hole bi n ary s tate described by 
McClintock fc RemillardI (|2006l ). The disc tem- 



date (see e.g. iBarnard et al.ll2011 
within) . 

3.3.2. XMM-Newton Rev 381 



The 0.3-10 keV EPIC-pn spectrum of RX J0042 
contained 21765 net source photons over ^^55 ks 
of live time, or ~0.4 count s^^. The detector was 
operated in Full Frame mode, with 73.4 ms frame 
time; hence, pile-up was negligible. 

An absorbed power law model failed to fit 
the spectrum; the best fit model yielded A^h 
~1.5xl02i atom cm'^ and F ~1.7, but xVdof 
= 412/348, with a null hypothesis probability of 
0.011. 

We also tried an absorbed disk blackbody 
model, since this is characteristic of the ther- 
mal high state identified in black hole binaries 
(e.g. iMcClintock fc RemillardI l2006l ). The best 
fit column density was a factor ~3 lower than 
the Galactic line-of-sight density (6.7x10^° atom 
cm~^), hence we fixed it to this value. This re- 
sulted in an inner disk temperature of 1.5 keV, 
but x^/dof = 765/349 and a null hypothesis prob- 
ability of 3x10^33 Hence, RXJ0042. 34^4115 was 
clearly not in the thermal high state. 

A disk blackbody -I- power law model de- 
scribed the spectrum very well, with A'h = 
2.4±0.06x 10^^ atom cm~^, inner disk tempera- 
ture kTin = 2.24±0.15 keV, F = 3.0±0.6, and 
X^/dof = 343/346. Figure [S] shows the un- 



perature is rather higher than usual, but consis- 
tent with the G alactic black hole binary system 
GRS1915-fl05; iMcClintock fc RemillardI (l2006h 
provide a sample spectral fit for GRS 1915-1-105 
with kTin = 2.19±0.04 keV and F = 3.46±0.02. 

Some authors have claimed that such a model 
is unphysical, because the power law compo- 
nent exceeds the thermal component at low en- 
ergies; they argue that the observed soft excess 
is an artifact of the two component model (see 
e.g. [Roberts et al . 2005: Gongalves fc S oria 20061 ). 
These arguments assume a compact corona that 
can only access photons from the inner disk. How- 
ever, there is substantial evidence for extended 
coronae in X-ray binaries at high luminosities; 
the ingress times of photo-electric absorption dips 
2_l_^-|^-|^t.in high inclination binaries indicate coronae with 



diameters of ^20,000-700,000 km (laiurch||2001 



Church fc Batucihska-Churchll2004 ). while broad- 



ened emission lines in Chandra observations of 
Cygnus X -2 suggest a hot, d ense corona of up to 
~10^ km (jSchulz et al.ll2009[ ). Such coronae would 
have access to the soft photons from the outer re- 
gions of the disc as well as the hot photons from 
the inner disc. 

Indeed, fitting the spectrum with a more phys- 
ically motivated model (diskbb + COMPTT in 
XSPEC) yielded the same values for iVn and kTin 
as the disk blackbody -I- power law model, along 
with a seed photon temperature of ^0.02 keV, an 
electron temperature of ~40 keV and an optical 
depth ~0.2; x^/dof = 343/344; the electron tem- 
perature was unconstrained, and good fits were 
obtained for electron temperatures of 100 keV and 
300 keV also, typical for the BH low state. Tying 
the seed photon temperature to kTin resulted in 
an unacceptable fit: x^/dof > 453/346 (null hy- 
pothesis probability < 1 x lO^"'). These results 
are entirely consistent with an extended, optically 
thin corona, and reject a compact corona that only 
sees the innermost region of the disk. 
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4. Discussion 

Although the original evidence for RX J0042. 3+4115 
being a black hole binary was contaminated by ar- 
tifacts in the XMM-Newton data reduction soft- 
ware, we now have new evidence that it is indeed 
a black hole candidate. 

Our constraint on the B magnitude (Mb > 
24.8 at the 3a level) is extremely useful for inter- 
preting the system. RXJ0042. 3-1-4115 is clearly 
not located in our galaxy, and there is no back- 
ground galaxy with B < 28; hence we locate 
RX J0042.3-K4115 in M31. The 3a upper limit to 
Mb = —0.4, meaning that a low mass donor is 
most likely. 

We have observed two distinct spectral states 
from this system, consistent with the low/hard 
and steep power law black hole states. The low 
state was observed at a 0.3-10 keV luminosity of 
2.08±0.08xl0^* erg s~^, while the steep power 
law state was observed at 2.41±0.05xl0'^^ erg 
s~^. Since transitions from the low state oc- 
cur at L <0.1 L khh in the 0.01-1000 k eV band 
for NS systems (jGladstone et al.l 120071 ). and in 
the 15-50 keV band for neutron star and black 



hole systems (jTang et al.l 120111 ). we suggest that 
RX J0042.3-h4115 was near 0.1 LEdd during Chan- 
dra observation 1575. Such a system would re- 
quire a ~20 M0 black hole; this is larger than 
for any Galactic black hole binary, but smaller 
than the dynamically confirmed black hole in 
ICIO X-1, which has a best mass estimate of 
32 .7+ 2. 6 M^, and a lower lim it of 23.1±2.1 Mq 
(|Silverman fc Filippenkoll2008l ). 

RXJ0042. 3-1-4115 has been persistently bright 
for the last ^12 years; this is quite unlike 
the transient behaviour of most Galactic black 
hole LMXBs. The only known persistently 
bright black hole LMXB is GRS 1915-1-105, which 
has remained bright sin ce its discovery in 1992 
([McClintock et all [2OO6I ). The X-ray behavior 
of GRS 1915-1-105 is unmatched by any Galac- 
tic black hole binary, and may be explained by 
the primary spinning in the same direction as 
the accretion disk at extreme speed s (>98% of 
the maximum iMcClintock et al.ll2006l ). Such pro- 
grade spinning allows the last stable orbit to be 
significantly closer to the black hole than for a 
non-spinning black hole, resulting in a higher disk 
luminosity and also a higher kTi„ (jZhang et al 



19971) . 

The inner disk temperature for RX J0042.3-K4115 
in the steep power law state (2.24±0.15 keV) 
was remarkably similar to that of GRS 1915-1-105 



(2.19 ±0.04 keV for an example spectrum. iMcClintock fc Remillard 
20061) . Hence, the persistent X-ray emission from 
RXJ0042. 3-1-4115 may also be due to extreme 
prograde spin. However, we note that 5 out of 
the 6 black hole candidates that we have asso- 
ciate d with M31 globular clusters are persistent 
also (Barnard etalj 120081; iBarnard fc Kolbll2009t 
jBarnard et al, 20111 ): such systems are consistent 
with theoretical prediction s for tidal capture o f 
main sequence donor stars (iKalogera et al. 20041 ). 
or ultra-compact sys tems with degenerate donors 



( Ivanova et al. 2010f) . Therefore, several mecha- 



nisms can promote persistently bright black hole 
binaries. 

Finally we note that compact corona models 
where the seed photon energy is tied to the inner 
disk temperature were all rejected by the XMM- 
Newton Rev 381 spectrum of RX J0042.3-H4115. 
However, free fitting of the seed photon energy 
yielded a good fit at 0.02 keV; this result gives 
strong support for an extended corona in the inter- 
mediate/steep power law state, which is able to ac- 
cess the cooler photons in the outer disc. We drew 
the same conclusions for the co nfirmed black h ole 
+ Wolf-Rayet binary ICIO X-1 (|Barnardll2O10l ). 
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Fig. 1.— A detail of the HST ACS/WFC image 
from observation j9udl7010; the exposure time 
was 4360 s, and the F435W (~B band) filter was 
used. The white ellipse represents the 3cr uncer- 
tainty in the position of RX J0042.3+4115. North 
is up, east is left. The brightest star consistent 
with the ellipse is indicated by a diamond; it has 
a Vega B band magnitude of 25.4±0.2. 
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Fig. 3.— EPIC-pn 0.3-10 keV lightcurve of 
RX J0042.3+4115 from the 60 ks Rev 381 observa- 
tion. The background lightcurve is shown in grey 
for comparison. We see that RX J0042. 3-1-4115 is 
significantly variable over short time-scales, vary- 
ing by a factor 2 in <10 ks. 
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Fig. 2. — Long-term, calibrated luminosity 
lightcurve of RX J0042.3-H4115 from -12 years of 
Chandra ACIS observations. The luminosity var- 
ied by a factor —6. Circles represent luminosities 
from freely fitted spectra, while crosses represent 
"faint" spectra. Uncertainties are quoted at the 
la level. 
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Fig. 4. — Unfolded spectrum for ACIS observa- 
tion 1575, multiplied by energy. It is well de- 
scribed by an absorbed power law with TVh = 
1.0±0.2xl02i atom cm^^ and F = 1.46±0.06; 
X^/dof = 211/204. Such a spectrum is charac- 
teristic of the low/hard state seen in neutron star 
and b lack hole binaries (|McClintock fc Remillard 
20061 ). 
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Fig. 5. — Unfolded spectrum for the 60 ks XMM- 
Newton observation Rev 381, multiplied by en- 
ergy. It is well described by a disk black- 
body (dashed) + power law (dotted) emission 
model, suffering line-of-sight absorption. A^h = 
2.4±0.6xl02i atom cm^^ kTin = 2.24±0.15 kcV, 
and r = 3.0±0.6; xVdof = 343/346. Such a 
spectrum is characteristic of the steep power law 
state seen in black hole binaries; the temperature 
is higher than is generally observed for Galac- 
tic black hole binaries, but is consistent with 
GRS 1915-^105 (|McClintock fc RemiUard,2006il . 
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